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Chapter 1 

Theoretical Background 



1.1 p Decay 

The history of beta decay began in 1896 with the discovery of radioactivity by 
Henri Becquerel. In the years 1899 and 1900 he identified f3 radiation as one com- 
ponent of radioactivity, and demonstrated that (3 rays are composed of electrons, 
comparing it with the properties of cathode rays. Unfortunately the rapidity with 
which (3 rays were identified after their discovery did not lead to an equally rapid 
interpretation of the (3 decay process. At that moment physicists knew that certain 
substances have the same chemical properties but different radioactive properties 
(isotopes). As experiments improved, the interpretations of the (3 decay phenomenon 
became more and more confusing. Chadwick in 1914, using a magnetic spectrom- 
eter, discovered the continuous spectrum of the /3 particles. Thus it was known 
that the electron spectra had monoenergetic lines and a continuous component, in 
contrast with alpha and gamma ray spectra, which were known to consist of mo- 
noenergetic lines only. Chadwick further demonstrated that the majority of events 
were part of the continuous component, the rest being monoenergetic electrons from 
electron conversion. The interpretation of the continuous electron spectrum was the 
main subject of the discussion. In 1914 Rutherford thought that the (3 electrons 
were all emitted from the nucleus with the same energy, but lost different fractions 
of this energy in collisions with the surrounding atoms, depending on the source 
thickness. But the main point was made by Lise Meitner in 1922. She realized 
that a quantized nucleus could not emit electrons of continuous energy. The known 
features of a and 7 spectra were correctly interpreted as due to transitions of nuclei 
from one quantum state to another. Thus the continuous electron spectrum was a 
unique feature of (3 decay. 

The measurement of the mass of different isotopes was determinant to conclude the 
existence of a neutral mass inside the nucleus, but the common thought was that this 
neutrality was due to the neutralization of the proton with an atomic electron which 
falls into the nucleus. Nevertheless, the uncertainty principle only allows electrons 
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inside the nucleus with an energy greater than ~ 100 MeVQ (the maximum energy 
of the electron continuous spectrum is around 10 MeV). And consequently the 
electrons could not exist " a priori" in the nucleus. 

In 1927 a crucial experiment was performed by Ellis and Wooster, in which they 
measured the total energy released in the disintegration of a 210 Bi source inside a 
calorimeter thick enough to stop all the emitted electrons. The endpoint of the elec- 
tron j3 spectrum was known to be E o =1.05 MeV, and the mean energy, E, of the j3 
electrons was known to be 390 keV. The calorimeter should have measured a total 
energy of 1.05 MeV if Rutherford's reasoning was correct. In fact they observed E 
= 344 ± 34 keV, which corresponded very well with the mean energy of the emitted 
electrons. The experiment was repeated in Berlin with an improved calorimeter by 
Meitner and Orthman in 1930 and the result was _E=337 ± 20 keV. These results 
were conclusive, Rutherford was wrong, but difficult to interpret at the time. 
In /3 decay not only was energy apparently not conserved, but the momentum and 
angular momentum were not conserved either. Pauli in 1931 proposed the idea of 
a very penetrating neutral particle of very small mass and spin 1/2 being emitted 
simultaneously with the electron. Thus, the problem of the conservation of energy, 
momentum and spin was solved. This proposal was made before Chadwick's dis- 
covery in 1932 of the neutral elementary particle, the neutron, and it was Fermi 
who named Pauli's particle as the neutrino. In this way the process of emission of 
electrons from the nucleus was explained and scientifically accepted. 
In 1934 Irene and Frederic Joliot-Curie observed the emission of positive electrons 
(later called positrons) in radioactive decay. And in 1938 Alvarez observed the elec- 
tron capture from the nearest electronic orbits. 

This brings us to our present knowledge of /3-decay which we can describe as, 

1. - (3~ decay: In this process a neutron (n) of the nucleus is converted into 

a proton (p) due to the week interaction, emitting an electron (e~) and an 
antineutrino (u). 

n p + e~ + v (1.1) 

2. - (3 + decay: In this process a proton (p) of the nucleus is converted into a 

neutron (n) due to the week interaction emitting a positron (e + ) and a neutrino 

(")■ 

p — > n + e + + v (1.2) 

3. - Electron Capture: In the electron capture process, also called inverse (3 + 

decay, an electron (e~) from the inner shell of the atom is captured by the 
nucleus, converting a proton (p) into a neutron (n) and emitting a neutrino 

p + e — > n + v (1.3) 

1 This rough calculation was made considering the size of the nucleus ~ 1 [fm] , but now we know 
that the size of the nucleus is ~ 6 [fm] , and considering this value the uncertainty principle allows 
electrons inside the nucleus with energies ~ 16 MeV 
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1.2 Fermi Theory of Beta Decay 

In 1934 Fermi developed a theory of the (3 decay process to include the neutrino, 
presumed to be masslessj as well as chargeless, dealing with the calculation of the 
transition probability of the process of decay This cannot be done starting from 
any other theory. A complete new force had to be introduced to explain the f3 
transition which converts a neutron into a proton (or viceversa) and at the same 
time produces an electron (positron) and a antineutrino (neutrino). Such a force 
was introduced by Fermi, using the analogy with the electromagnetism. 
The transition probability in the /3 decay process can be given in terms of the first- 
order, time-dependent perturbation theory, later called the Fermi Golden Rule, 

A = T I(/I^K)| 2 P/ (1-4) 

where (f\Hp\i) is the matrix element of the beta interaction Hp between the initial 
state |i) and the final state |/) of the complete system (nucleus and other relevant 
light particles), and pf is the density of final states in the final system. The final 
state of the system is specified by the electron and neutrino momenta and energies, 
(p e ,E e ) and (p u ,E u ), with E u = cp u . If Ej denotes the energy in the final system, 
we have 

2 p 2 e dp e dVL e pldp u dVL u 
PfdEf = V (2^)3 (2^)3 (L5) 
where V is a spherical volume in the momentum space and 

E f = E e + E v = E e + p u c (1.6) 

Without fixing the directions of the momenta, we have 

" = *Wf^ (1.7) 

Fermi did not know the mathematical form of the interaction Hp in /3 decay. But in 
considering all possible forms consistent with special relativity, he showed that Hp 
could be replaced by an operator O x which could, mathematically, take the form of 
a vector(V), scalar (S'), pseudoscalar(P), axial vector (A) or tensor (T). So, for beta 
decay, 

(f\H P \i) = w fi = E^/ d ^M^:] o x ^ 



where $/, ip e and ip u are the final wave functions of the nucleus after the decay, the 
electron and the neutrino. The value of g x determines the strength of the interaction. 
The electron and neutrino are treated like free-particles, thus their wave functions 
have the form, 

V> e (r) = -j==e ip ~^ n , Mr) = -j=e ip ^ /n (1.9) 



2 Today we think that the neutrino has some mass, although very small 
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If we expand the exponentials, using the fact that over the nuclear volume pr -C 1, 
we have the allowed approximation. If we replace the electron and neutrino wave 
functions in eq jl.81 and use the allowed approximation, the matrix element is now 
gMfi = g f dr§*fO&i, so the decay rate is 

g 2 f dp v 2 2 



A = »^ |M '''*« (L10) 

for a fixed E e the term = 1/c. On the other hand, if we define Q as the total 

decay energy, the momentum of the neutrino is p v = (Q — T e )/c. The decay rate is 
now 

A = / d P \M fi \Y(Q - T e f (1.11) 

but here we are not taking into account the interaction between the beta particle 
and the Coulomb field in the daugther nucleus. In order to take this effect into 
account an additional factor was added to the decay rate, F(Z',p) where Z' is the 
atomic number of the daughter nucleus. The total decay rate is now 

X = tWt? / - T e ) 2 F(Z',p)dp. (1.12) 

2ir 6 n c 6 Jo 

This integral only depends on Z' and the maximum electron energy E , and it is 
represented as 

1 rPmax 

f(Z',E ) = — / p\Q-T e fF{Z',p)dp (1.13) 
mid J 

where the constants have been included to make / dimensionless. This function is 
known as the Fermi integral and it is tabulated for values of Z' and Eq. 
As A = ln(2)/Ty 2 , we have, 



g 2 m 5 e c 4 \M fi 



mn = .r_»T.; „ (i-") 



The equation 11.141 gaves the comparative halflife or ft value. 



1.3 Fermi and Gamow Teller Decay 

In the allowed approximation the electron and the neutrino are created at the 
origin (r = 0). As the orbital angular momentum is zero(/ = 0), the only change in 
the angular momentum is due to the spins of the electron (s e = 1/2) and neutrino 
(s v = 1/2). If the two spins are antiparallel (S = 0), there is no change in the 
nuclear spin: AI = |Jj — = 0. This is known as Fermi Decay (F). If the electron 
and neutrino spins are parallel {S = 1), they carry away an angular momentum of 1, 
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therefore F and If must be coupled through a vector of length 1: 7j = If + 1. This is 
possible for AI = or 1 (except for F = and 7/ = 0, when only Fermi transitions 
can contribute). This is known as Gamow-Teller Decay {GT). They are governed 
by different operators which we will called Op and Oct m an obvious notation. 

The matrix elements of F and GT decays must be written separately. The matrix 
elements for Fermi Decay, following equation ll.8[ is 

M fi = M F = g v J dWfOj&i (1.15) 

and for the Gamow-Teller Decay 

"» = ** = "j*«P<**< (i,6, 

The Fermi operator is expressed as Op = J2t=i T ±(^) where t + (t_) is the 
isospin ladder operator converting the proton (neutron) wavefunction into a neu- 
tron(proton) wavefunction. For Gamow-Teller decay the corresponding operator is 
Oct = Y^t=i ^(*) r ±(*) where a(i) are the Pauli matrices, which also act on the ith 
nucleon. Now we can rewrite the equation 11.141 in terms of the Fermi and Gamow- 
Teller matrix elements. 

K K 2ln(2)n 3 h 7 

jl 1/2 = 2 = 2 ) where K = — 5 — — — 

|M F | 2 + %|M GT | 2 B(F) + 9 -4-B(GT) 9 v m l c 

9v 9v 

(1.17) 

and B(F) and B(GT) are the dimensionless Fermi and Gamow-Teller strengths. 



1.4 Isotopic Spin of the Nucleon 

The charge symmetry and charge Independence of the nuclear interaction and 
the almost equality of the masses of the neutron and the proton strongly suggest 
that they are the same particle (the nucleon) in two different charge states. [U El E] 
In analogy with the reorientation of the two degenerate states of the electron spin 
under the effect of a magnetic field in the real space, a nucleon under the influence 
of a electromagnetic field, can differentiate the two degenerate isospin states or 
charge states in the charge space or isospin space. A nucleon with third isospin 
component T z = —1/2 (down) is defined as a proton, and a nucleon with third 
isospin component T z = +1/2 (up) is defined as a neutronJl. The isospin obeys the 
usual rules for angular momentum vectors. The isospin vector r has a length of 

3 The isospin assignment for a neutron and a proton, which is adopted here, is conventional in 
nuclear physics and has the advantage that the heavy nuclei with large neutron excess have their 
isospin aligned in the direction of the positive z axis in the isospin space. In elementary particle 
physics, the inverse sign convertion is usually adopted. [3 
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t(t + Vj% and with 3-axis projections r 2 = m T h. 
For a system of several nucleons, the isospin follows the same coupling rules as the 
ordinary angular momentum vectors. The 3-axis component of the total isospin 
vector, T z , is the sum of the 3-axis components of the individual nucleons, and for 
any nucleus 



A 1 

T z = Y,t^) = -(N-Z) (1.18) 



i=l 



expressed in units of ft. [5] 



1.5 Isospin Symmetry 

The assumption that the attractive nuclear force is independent of the charge of 
the interacting nucleons is one of the bases of our understanding of nuclear physics. 
This assumption hides a simple and elegant symmetry: The isospin symmetry (or 
charge symmetry). 
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Figure 1.1: Schematic illustration of the isospin symmetry transitions in the A = 54 
system. The level schemes, neglecting Coulomb displacement energies, are shown to 
illustrate the isospin symmetry structure of states in different T z nuclei. 

In the isospin approach, all nucleons must have an isospin t = 1/2 differentiating 
protons and neutrons by its projections. The use of isospin gives us a simple way to 
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classify nuclei according to the third component of their isospin. 
A state in a nucleus with N = Z corresponds to a particular configuration of neu- 
trons and protons. If an identical state cannot be constructed in the neighbouring 
nucleus with N — 1 neutrons and Z + 1 protons by exchanging a neutron by a proton 
then the original state must have isospin T = 0, because it can only be constructed 
in a nucleus with N — Z = 0, which has projection T z = 0. The N = Z nucleus 
can also have a state with T = 1, which, because T z ± 1 is also allowed, can be 
constructed identically in the N ±1, Z =p 1 systems(see Figure [TTTj) . [I] 



1.6 Experimental Studies of Isospin Symmetry 

The T z = +1 to T z = transitions could be studied in charge exchange (CE) 
reactions ((p, n)-type [6], [7] or ( 3 He, t)-type[8],[9]), and the T z — — 1 to T z — tran- 
sitions can be investigated in (3 decay experiments as shown in Figure 11.11 Such 
studies have only been performed in detail so far for light nuclei, such as the A = 26 
system ( 26 Mg, 26 Al and 26 Si)[TU]. One major constraint for these studies in light 
nuclei is that (3 decays can access only a few low-lying states due to the restriction 
imposed by the Q/3-window. Among the T = —1 — > candidates, the analogue 
transitions in the fp shell nuclei are well suited to an accurate study of isospin sym- 
metry. The reason is that due to relatively large Qec values, the (3 decay studies 
should allow B(GT) measurements up to high excitation energies in the daughter 
nuclei. 



1.6.1 Charge Exchange reactions (T z = +1 — > T z = 0) 

It is known that the CE reactions, such as (p, n) or ( 3 He,£) reactions, at inter- 
mediate incident energies are valuable tools for B(GT) studies up to high excitation 
energies in the final nucleus, since they are not limited by the Q window as in /3 
decay. It has been found that CE reactions performed at angles around 0° and 
intermediate energies (E in > 100 MeV/nucleon) are good probes of GT transition 
strength due to the approximate proportionality between the cross section at 0° and 
the B(GT) values [II], 0, 



where J aT (0) is the volume integral of the effective interaction V ar at momentum 
transfer q = 0, K is a kinematic factor, N aT is a distortion factor, and <Jgt(0°) is 
the unit cross section for a GT transition at 0°. CE reactions have greatly improved 
in recent years, because the poor energy resolution in the pioneering (p, n) work has 
now been overcome by the use of the ( 3 He,t) reaction [T2"].[T5]. 



_d_ 



a CE (0 o ) ~ KN aT \ J aT (0)\ 2 B(GT) = a GT (0°)B(GT) 



(1.19) 
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1.6.2 p Decay (T z = -1 T z = 0) 

The transition strengths B(GT) are directly obtained from /3 decay ft values 
which are derived from measurements of the Qec value, the half-life and the branch- 
ing ratio of the transition of interest. 

In a (3 decay, the partial halflife ti/ 2 multiplied by the /-factor is related to the 
B(GT) and the Fermi transition strength B(F), as we saw in eq. 11.171 which has to 
be corrected by a Coulomb factor, 

/Tl/2 = ^5(F)(1-5 C ) + A^(GT) (L20) 

where i^=6144.4±1.6, A=(7A/5 , y = -l-266±0.004, 5 C is the Coulomb correction factor. 
The Fermi strength can be calculated theoretically, and its value is B(F) = \N — Z\, 
on the contrary, the B(GT) has to be taken from experiment. 
Uncertainties in the experimental B(GT) values originate from uncertainties in the 
decay Q-value, the total halflife Ti/ 2 and the branching ratios, determining t. The 
accurate determination of the feeding ratios to higher excited states is more difficult 
due to the smaller /-factors. On the other hand, in the ( 3 He, t) measurements 
studying the corresponding GT transition, relative transitions strength to these 
higher excited states can be obtained accurately. The feeding ratios in the (3 decay 
can be deduced using these values and /-factors that are calculated using the decay 
Q-value. These feeding ratios can then be converted into absolute B(GT) values by 
the normalization process using the total halflife of the f3 decay. 
The inverse of the total /3-decay halflife T1/2 is the sum of the inverse of the partial 
halflife tp of the Fermi transition to the Isobaric Analogue State (IAS) and the 
partial halffifes (tj) of GT transitions to the zth GT states. 

1 1 1 , 

i/z r i=GT 1 

Applying eq. I1.20l one can eliminate both tp and the t,, so that 

^- = )? [B(F)(l - 5 c )f F + * 2 MGT)fA (1.22) 

Tl /2 K V i=GT J 

where fp and /j are the /-factors of the /? decay to the IAS and to the ith GT 
states, respectively, Bi(GT) is the B(GT) value of the transition to the ith GT state, 
and B(F) = \N - Z\. 

In order to relate the strengths of GT and Fermi transitions in a CE reaction, the 
ratio R 2 of unit GT and Fermi cross sections at 0° it is introduced. 
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Due to the isospin symmetry, this ratio R 2 is expected to be the same for all the 
T z = ±1 — > transitions. Eliminating Bi(GT) by using R 2 , we get 



B(F)(1-5 C ) 



1/2 




(1.24) 



From this merging method point of view, B(F) and 5 C depend on which nuclei 
we are working with, a F and a GT are obtained from CE reactions, fa and fp are 
obtained from calculations and the total halflife Ti/ 2 is obtained from (5 decay. 

In summary, the relative values of Bi(GT) can be obtained from the CE reaction. 
Consequently, if we have an accurate value for T ± / 2 from (3 decay studies, we can use 
the CE information to extract absolute B^GT) values to all the states observed in 
the reaction. 
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Chapter 2 

Experimental Setup 



The experiment to measure the 54 Ni beta decay was carried out at the Leuven 
Isotope Separator On Line (LISOL) at Centre de Recherche du Cycloron at Louvain 
La Neuve (Belgium) in April 2006. 



2.1 The Production 

The 54 Fe target located in a gas cell was bombarded with a primary 3 He beam 
at 45 MeV coming from the CYCLONE110 cyclotron. Several radioactive nuclei 
were created in this reaction but since we have an on-line mass separator we are 
only concerned about the production, selection and mass separation of the 54 Ni and 
54 Co nuclei. 

The decay of 54 Ni is difficult to study because its production is always over- 
whelmed by the direct production of its daugther nucleus 54 Co. Indeed the calcu- 
lated production cross section for the 54 Ni has a maximum of 30 /zbarn at an energy 
of 45 MeV for the 3 He beam. At the same energy the 54 Co has a production cross 
section of 5 mbarn (166 times bigger) [TP] • 



2.2 The LISOL Separator 

The 54 Fe( 3 He,2np) 54 Co and 54 Fe( 3 He,3n) 54 Ni reactions are produced in the gas 
cell. A laser beam selectively ionizes either Ni or Co atoms. The resonant ions are 
extracted and accelerated by a DC electrical field and directed to the mass separator 
which separates the single- charged, A=54 mass ions by A/q. In the following we 
give a detailed explanation of the different components of the LISOL device. 
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Cyclotron 




Figure 2.1: Leuven Ion Separator On Line. 

2.2.1 Gas Cell 

The gas cell is a gas catcher source and is the first component of the so-called 
Laser Ion Source (LIS). The nuclear reaction occurs in the gas cell filled with a 
noble gas flowing through it. The reaction products are stopped, thermalized and 
neutralized in noble gas. This avoids chemical reactions with the nuclei produced. 
The body of the cell is made of stainless steel and is electro-polished to reduce the 
level of roughness of the surface. 

In order to optimize the efficiency and the selectivity of the ion of interest, a LIS 
has to fulfill some important requirements such as (a) The reaction products, recoil- 
ing out the target, have to be stopped in the buffer gas, (b) all reaction products 
must be neutralized before arriving in the laser ionization zone, (c) avoid as much 
as possible the interaction of the ion of interest with impurities in the buffer gas to 
avoid the formation of neutral molecules, (d) the evacuation time should be shorter 
than the mean diffusion time to the wall in order to avoid sticking of the reaction 
products to the wall of the gas cell, (e) the evacuation time should be shorter than 
the lifetime of the isotope of interest, (f) the time between two subsequent laser 
pulses should be shorter than the evacuation time of the ions produced in the laser 
ionization region and (g) the survival time of the laser produced ions has to be longer 
than their evacuation time. 



2.2 The LISOL Separator 
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2.2.2 Resonant Laser Ionization 

After the thermalization and neutralization of the reaction products we need to 
select and extract the ion of interest. The resonant laser ionization technique is used 
to ionize a specific element. 

Two tunable dye lasers are pumped by two excimer Xe-Cl (308 nm) lasers. The 
atoms thermalized in the ground state level are excited in two steps. The first step 
laser Ai excites the atom to an intermediate level. The second step laser A2 may 
ionize and excite atom to the continuum, to an autoionizing state or to a Rydberg 
state and where it is ionized by collision with the buffer gas (See Figur d2~2"j) . 
This two dye laser beams are focused by lenses and directed by prisms to the gas 



Autoionizing States 



A*+e 



First Ionization 




Figure 2.2: Atomic level scheme and the different ways to ionize an atom. 

cell located at 15 m from the optical setup. In order to focus the two beams at an 
optical distance of 15 m a screen is located at that distance but not in the same 
position as the gas cell. 

To choose the most efficient resonance wavelength of the two dye laser beams, a 
fraction of each beam is deflected and directed to a reference cell in a vacuum cham- 
ber. An atomic beam is produced by evaporating the element investigated from a 
resistively- heated crucible. 

After the laser ionization, the ionized atoms and compound molecules [1 5 j are at- 
tracted by a potential difference and ejected from the gas cell. 
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2.2.3 Sextupole Ion Guide 

Once the selected atom is ionized it can re-combine with the impurities (nor- 
mally H 2 0) in the gas cell forming molecules [TH]. These compound molecules are 
heavier than a single atom so the isotope of interest is lost in the subsequent mass 
separation step. This problem can be solved if a sextupole ion guide is located at 
the exit hole of the gas cell [16] . 
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Figure 2.3: Sextupole Ion Guide (SPIG) 



The SextuPole Ion Guide SPIG, as shown in Fig.f l2.3p consists of six silver 
rods (1.5 mm diameter and 124 mm lenght) cylindrically mounted and forming a 
hexagonal structure parallel to the beam path. An oscillating radio frequency voltage 
(V r /) is applied to the rods, with every rod in antiphase with the neighbouring rods. 
The amplitude and the frequency of the rf signal are chosen in order to confine 
the ions in the centre of the SPIG. Hence the gas coming out of the gas cell can 
expand through the gaps between the rods, while the ions are radially confined and 
transported to an extraction electrode with the velocity of the gas jet. 
In addition to the rf voltage, a DC voltage is applied between the rods and the 
gas cell exit (grounded). Hence a molecule could dissociate and liberate the ion of 
interest, thus increasing the selection production efficiency of the LIS. 

After that, the selected ion is mass-separated in a mass separator. The mass 
separator is a 55° dipole magnet of 1.5 m radius. The ions are separated according 
to the mass-over-charge ratio since they are, in general, in a single charge state, 
normally 1 + . Hence they are separated in mass. The ion beam of the mass selected 
is implanted into a tape located in the middle of two Miniball germanium detectors 
(see Figure P2T31 . as described in the next chapter. 



2.3 Radiation Detectors 
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2.3 Radiation Detectors 

In the (3 decay of 54 Ni we expect to have three type of radiation: betas, neutrinos 
and gammas. It is difficult to detect the neutrinos but we can make a sensitive setup 
to detect radiation and 7-rays. For the 7-ray detection we used two MINIBALL 
Ge cluster detectors and for the j3 particles we used plastic detectors. Three 1.3 mm 
thick plastic (3 detectors were placed around the implantation point (tape ending, 
see Figure 12731) . In most of the cases the (3 particle only deposits part of its energy 
in the thin plastic. These detectors were used to determine if a (3 particle passed 
through the detector or not. No energy was registered. 




Figure 2.4: Top and side view of the setup of (3 and 7 detectors at Louvain la Neuve 
facilities. 

Two MINIBALL clusters were located as close as possible to two of the scintilla- 
tors. The MINIBALL detector is a cluster of three hyper-pure germanium crystals, 
each of them six-fold segmented. The crystals and the capsules are tapered in a 
hexagonal shape to fit three of them side- by-side in a single cryostat. A positive 
high voltage is applied to a lithium-diffused centre contact, which penetrates deep 
inside the germanium. Six boron-implanted contacts on the surface of each crys- 
tal constitute the segmented outer contacts. An incident 7-ray interacts with the 
Ge crystal and transfers its energy either totally (photoelectric effect) or partially 
(Compton or pair creation) to an electron, this electron loses its kinetic energy, 
creating a number of electron-hole pairs. Within the fully depleted region of germa- 
nium between the contacts, the charge created goes to the corresponding electrode. 
Thus, the electrons go to the centre contact (core signal), and the holes go to the 
outer contacts (segment signal) (See Figure [2~B]) . 
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The charge collected on the crystal contacts is integrated by charge- sensitive resistor 
feedback preamplifiers, developed and produced by the Institute of Nuclear Physics 
in Cologne, a member institute of the MINIBALL Collaboration. [T8] 



Miniball Cluster 



Cores 




Figure 2.5: Schematics showing cores and segments in a Miniball Cluster 



2.4 The Tape System 

After the mass separation, the ions are implanted on a thin 1.25 cm wide 
aluminium-coated mylar tape. The thickness of the tape is chosen to stop the 
beam completely. The tape rewinding system is located below the detectors. The 
tape from one of the two coils goes up, intercepts the beam path, goes up again 
and turning fully down, comes again below the detector level and finally ends in the 
second coil. The part of the tape where the ions are implanted is moved down to 
a safe distance away from the detectors, to avoid the background produced by the 
activity of long-lived ions close to stability. The whole tape system is kept in the 
same vacuum as the beam transport line. There are three thin mylar windows in the 
chamber around the implantation area to reduce the attenuation of the radiation of 
interest. 

A schematic picture of the detector setup and the tape system is shown in Figure 

EM 



2.5 The Data Acquisition System 
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Figure 2.6: Schematic view of the detector and type system setup 



2.5 The Data Acquisition System 



The preamplified signals are provided to the analogue inputs of the DGF-4C 
modules. The DGF (Digital Gamma Finder) can measure the amplitude and the 
shape of the incoming pulses simultaneously. The energy information is extracted 
directly from the preamplified signal height and the time information is obtained 
from a built-in high-precision internal clock after the event validation. 
In our detector system we have 21 signals from each MINIBALL 7 detector (3 sig- 
nals from the cores and 18 segments) and another 3 signals from the three plastic 
scintillator (3 detectors. Since an incident 7 ray is registered by the core and some 
of the segments of the same crystal, there is a time correlation between the signals 
from one crystal. Further on there should be no segment signal unless there is a 
signal from the corresponding crystal core. For these reasons, the DGF modules are 
split into trigger groups, in which each group represents one crystal and the core 
channel is the main trigger of the group. The logic signals from the (3 detectors are 
provided to one DFG module and their inputs are set independently. To estimate 
the total acquisition live time, a pulser signal feeds one channel of a DGF module. 
This value is represented by the ratio between the registered pulser signals to the 
number of pulses provided by the pulser. 
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2.6 Data Analysis 

A new C++ software code has been created by Oleg Ivanov and Dieter Pauwels 
from Katholieke Universiteit Leuven, to read and analyze the data acquired by 
the detection setup. The object oriented ROOT package developed at CERN is 
used for building the output histograms. While building the singles spectra is a 
simple procedure, the ft-nj and wy coincidences, where n represents the number of 
7-detectors (crystals), are based on the event time stamps. These coincidences are 
stored in ROOT trees for further analysis [18J. 



Chapter 3 

The Experiment 



The aim of this experiment was to study the GT transitions in the /3-decay of 
54 Ni analogous to the GT transitions of 54 Fe( 3 He,t). The main objectives are to 
measure 1) the total halflife Ti/ 2 of 54 Ni with high accuracy, 2) the decay branching 
ratios of the 54 Ni to the ground state and to the first GT state and, if possible 3) 
the decay to at least one extra GT excited state in 54 Co. 

Figure [37T1 shows the known decay scheme of 54 Ni and 54 Co 9,m . The /3 + decay of 
54 Ni feeds the ground state of 54 Co and at least one exited state (I n =l + ) with an 
energy of 937 keV. The decay of the 54 Co m is accompanied by a triple 7 cascade of 
411 keV, 1130 keV and 1408 keV 7 rays. The 54 Co ground state decays 100% to the 
ground state of 54 Fe. 



3.1 Calibrations 

The resolution check for the MINIBALL detectors was performed with a 60 Co 
source. The measured resolutions for all MINIBALL crystals at 1332 keV are shown 
in the next table for the "Heidelberg" and "Leuven" clusters. 



Heidl.A 


3.4 keV 


LeuvenA 


3.3 keV 


Heidi. B 


2.7 keV 


LeuvenB 


3.1 keV 


Heidi. C 


2.8 keV 


LeuvenC 


3.0 keV 



The 7 efficiency calibrations (See Figure I3.2p were performed after the experi- 
ment with a 152 Eu positioned on the tape at the point of implantation. The activity 
of the source at the time of the experiment was 23.87 kBq. 

The same 152 Eu source was used before the experiment, for the on-line analysis 
energy calibration. Afterwards, for off-line analysis, an internal energy calibration 
was performed using the 54 Co m gamma cascade present in our experiment for each 
ROOT file created after more or less each one hour of running time. 
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Q= 8800(50) keV 



Figure 3.1: 54 Ni and Co 9 ' 771 decay scheme. 

3.2 Control Measurements during the Experi- 
ment 

With all the calibrated spectra, several control measurements were performed 
before and during the experiment in order to optimize the production and mass 
separation of 54 Ni and 54 Co s ' m . 

1. - Laser tuning on Co: To optimize the selectivity of Co ions by laser resonance 

in the gas cell, we looked at the ratio in intensity of 411 keV and 511 keV 7-ray 
lines from the decay of 54 Co 9 ' m in the singles and — 7 coincidence spectra 
with a fixed mass 54 in the mass separator. The cyclotron current and the tape 
steps were also fixed. Short measurements of 100 seconds with and without 
laser resonance were made. 

2. - Mass tuning: Once the selectivity was optimized, we proceeded to tune the 

mass selection on the separator with laser resonance on Co looking at the 
integral of the 411 keV line in the singles and — 7 spectra, taking short 
measurements of 100 seconds. 

3. - Gas Buffer Pressure and Cyclotron Current: With the mass optimized 

and with laser resonance on Co, the optimal gas buffer pressure was found 



3.2 Control Measurements during the Experiment 
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Figure 3.2: 7 Efficiency measurement of the MINIBALL detectors with a 152Eu 
source after the experiment. 



looking again at the 411 keV line in the singles and f3 — 7 spectra, taking short 
measurements of 100 seconds. After this the cyclotron current was tuned 
following the same procedure as described above. 

4.- Tape Step and Implantation-Decay Macrocycle: In the previous 54 Ni 
f3 decay experiment developed by I.Reusen et al.[T7],[T2] the measured halflife 
was 106(12) ms. Consequently, a proper 400 ms and 600 ms implantation- 
decay macrocycle was chosen for this experiment. The tape was moved one 
step every six macrocycles in order to minimize the time lost during the tape 
rewinding. 

In the following table the starting settings are summarized. 



Mass 


53.94 


Cyclotron Current 


4.0 fiA 


Gas Buffer Pressure 


760 mbar 


ID Macrocycle 


400 ms / 600 ms 


Tape Movement Period 


6 Macrocycles 



These settings were checked and, if necessary, modified during the experiment in 
order to increase the selectivity and mass separation of the ion of interest. 
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3.3 Data Analysis 

After seven days of experiment, the data acquisition time with laser resonance 
on Ni was 115 hours. The time used for the laser resonance on Co was 16 hours. 
The total number of counts in the 411 keV 54 Co m decay radiation in the laser on 
Ni total spectrum was 94,817 while the number of counts in the 411 keV peak in 
the laser on Co total spectrum was 46,903. Finally, the settings adjustment time 
was 2.5 hours. Singles and /3-7 coincidence spectra were constructed in ROOT trees 
to facilitate the on line analysis. These pre-processed data were also used for the 
off-line analysis. 

In general the data analysis was done with the /3-7 coincidence spectra. As we are 
interested also in the high energy GT excited states, a veto condition (see Figure 
13. 3j) was imposed in order to reduce the background produced by the penetration of 
the /3 particle into the MINIBALL crystal. The main idea is that if a (3 particle is 
detected in a plastic detector, the gamma detector situated behind it is not used. In 
Figure 13.41 we see the difference between the spectrum with the no- veto condition, 
which has 5,984 counts in the 937 keV peak and with the veto condition, which 
has 3,886 counts in the 937 keV peak. Although the upper spectrum has higher 
statistics, the lower one has a better peak-to-total ratio. 

Figures 13.51 and 13.61 represent the total 7 spectra of mass 54 activity with the 
laser set on resonace for Ni and Co respectively. Both of them were constructed with 
the no/3 penetration veto condition. In the laser on Ni spectrum we can observe 
the 937 keV 7 line known in the decay of 54 Ni. The gamma cascade lines 411 keV, 
1130 keV and 1408 keV from the 54 Co m decay and the 922 keV sum peak (411 keV 
+ 511 keV) are also present. Another peak is observed in the spectra at 683 keV 
with a measured halflife of ~ 180 ms, but it is not yet identified. See also Figure I3~7H 
for the low energy part of the spectra. In the laser set on Co spectrum, the same 
cascade lines and sum peaks are present but not the 937 keV line. 



3.4 Halflife Analysis 
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3.4 Halflife Analysis 

From the previous experiment [17] we know that the 54 Ni halflife is 106(12) ms, 
thus if we take into account the macrobeam structure (400 ms of implantation and 
600 ms of decay) and the microbeam structure due to the mass separator time 
behaviour during the implantation, we can model our Implantation Decay Curve, 
simply by solving the decay differential equations. 



^(t) = P(t) - X 1 N 1 (t) N 1 (0) = (3.1) 
j t N 2 (t) = X 1 N 1 {t) - X 2 N 2 (t) N 2 (0) = (3.2) 

4 

P{t) = a(H{t - 50(2i - 1)) - H(t - 100i)) 

i=l 

where N\ is the 54 Ni activity, N 2 is the 54 Co activity, Ai and A2 are the respective 
decay constants and P(t) is the implantation rate of 54 Ni, taking into account the 
microbeam structure. H(t) is the Heaviside step function and a is the constant 
implantation rate when P(t) ^ 0. 

Solving the differential equation with the initial conditions for N\{£) and N 2 (t) using 
Maplel0' m , we obtain the implantation decay curves of 54 Ni and 54 Co (Fig J3.8l) . 



On the other hand, from our experiment we obtain a f3 gated 7 spectrum, thus 
we can reproduce the ID curve looking at the time evolution of the 937 keV peak 
during one macrocycle. The time difference between each of the spectra (50 ms) 
was chosen in order to be consistent with the time behaviour of the mass separator. 
This way we can find the halflife value of 54 Ni by looking at the time behaviour of 
the integral of the 937 keV 7 line. 

The analysis of the halflife of 54 Ni can be deduced from the time behaviour of 
the area of the 937 keV peak after background substraction in the laser on Ni 7 
spectrum. 

As a first approach to obtaining the halflife we integrated a fitted Gaussian with 
937 keV centroid, subtracting a linear background close to the Gaussian (centroid ± 
twice FWHM region). Figure 13791 shows the evolution of the 937 keV peak during an 
ID macrocycle. The last four spectra corresponding to the last 200 ms show that the 
intensity of the 937 keV line became of the same order as the background. For that 
reason, the last four points were neglected in the fit of the micro implantation-decay 
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curve (see Figure l3.10p . In this case, the error obtained in fitting a Gaussian for each 
experimental point in Fig J3.10l is calculated with the ROOT minimization package 
MINUIT.0 The value obtained for the halflife of the 54 Ni in this case was Ti/ 2 = 114.9 
± 6.0 ms. The halflife error was obtained from the fit to the implantation-decay 
curve including the micro structure, again using the ROOT minimization package 
MINUIT. 

Another attempt to deduce the halflife as precisely as possible was made by fit- 
ting two Gaussians, one on the 922 keV sum peak and the other on the 937 peak, the 
second one with a fixed width (see Figure 15. lip . The value of the width was taken 
from the width of the Gaussian fit in the total spectra. The so called long back- 
ground (linear background taken from 800 keV to 1100 keV) was used to subtract 
the background in the region of interest, in order to reduce background fluctuations 
in the vicinity of the 937 peak. A reduced error bar is obtained with this method 
in comparison with the first attempt. Again, the error obtained by fitting a dou- 
ble gaussian for each experimental point in Fig J3.10l is calculated with the ROOT 
minimization package MINUIT. However it is also clear from this fit that the last 
points do not behave in the expected way and should therefore be ignored in the Ti/ 2 
fit. This tells us that due to the lack of statistics (peak indistinguishable from the 
background) at the end of the macrocycle, we can eliminate those points in order to 
fit a proper curve to determine the halflife. The value of the halflife of the 54 Ni for 
this case was T\/2= 113.8 ± 4.4 ms. The halflife error was obtained from the fit of 
the implantation-decay curve including the micro structure using again the ROOT 
minimization package MINUIT. 

A third method to measure the halflife of the 54 Ni was to integrate directly the 
937 keV peak and subtract the long linear background (See Figure l3~T3]) . The er- 
ror bar is a statistical error obtained from the integration of the 937 keV peak. If 
we again neglect the last 200 ms of the macrocycle, the halflife value was T\/2= 
113.7 ± 4.6 ms(See Figure I3.14p . The halflife error was obtained from the fit of 
the implantation-decay curve including the micro structure using again the ROOT 
minimization package MINUIT. 

In spite of the good agreement shown between these three different halflife anal- 
yses, it is worth mentioning that the estrange behaviour at the end of the ID curve 
clearly suggested that some correlation existed between the range of the background 
chosen and the halflife obtained. Indeed in the 937 keV region, in addition to the 
922 keV sum peak, there is the Compton edge of the 1130 keV 7 peak. This made it 
very difficult to decide if the background that we were subtracting was the correct 
one or not. 



1 MINUIT calculates an error matrix with the parameters. This error matrix is the inverse of the 
second derivative matrix of the best parameter values (the function minimum) and the diagonal 
elements of the error matrix are the square roots of the individual parameter errors, including the 
effect of correlation with the other parameters. 



3.5 An Attempt to Reproduce the Shape of the Background 
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3.5 An Attempt to Reproduce the Shape of the 
Background 

A semi-empirical model for the 7-ray response of the germanium detectors [2U] 
was used to estimate the correct shape of the background for background subtraction 
under the 937 keV line. In this analysis was include: 

a.- The full energy peak: A Gaussian distribution with variance B\ and nor- 
malization parameter A\ 



b.- The Flat continuum: A constant function when folded with a Gaussian 
function gives an expression of the form 



c- The single Compton scattering: The expression for the energy deposited 
in a single Compton scattering is well documented in textbooks. In order to 
account for the roll-off at the Compton edge that is observed in experimental 
spectra, Y.Jin et al. determine the electron energy function with a comple- 
mentary error function, 




(3.3) 




(3.4) 



f 3 (E ,E) = A 3 
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f 3 (E ,E) = 



E> E, 
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(3.5) 




and 



E 



1 + 



mpc 2 
2E 



d - Multiple Compton scattering Y.Jin et al. have derived expressions for 
the energy deposition due to two and three Compton scatters. The derivation 
assumed that the spatial distribution of once-scattered and twice-scattered 
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photons within the detector, are uniform. They give a set of equations for the 
energy deposition function for the double Compton scattering f 4 {E G , E) 



f 4 (E ,E)=0, a">a 

POL „ 

/ Ol 

U(E ,E) = A 4 H(a,a')H(a',a")da', - 

J a" I 

f 4 (E , E)=A 4 / 1 ~ 2a " H{a, a')H{a', a")da', - 

J - °L 1 



< a" < a 

2a 



a „ a 



< a < 



where 



+ 4a l + 2a 

a 

U(E ,E) = 0, 0<a"< TT ^ (3.6) 



H(a, «) = —+ - + — -— - + ~n2 + — + — 

cr \a cr a 6 ) \a ) a z a A a 6 cr 

and 

a — -, a = -, a = -, E = h -E. 

m c z m c A m c z 

The energy deposition function f^Eo, E) for the triple Compton scattering is 
given by the following set of equation 



f 5 (E ,E) = 0, a">a, 



f 5 (E , E) = A 5 I" [ H(a, P)H(P, 7 )ff ( 7 , a")<MA - 

J a" J a" l 



f 5 (E ,E) = A 5 j a J /f(c^)ff(A 7 )ff (7,^7^ + 



a 



+ 2a 



< a < a 



l+2a 
^ ft 



a 



A 5 I r_2" H(a, p)H(0, 7 )/f ( 7 , a»)d^ Th~4« < °" < TT~2a 

l-2a" 1+2/3 



tt tt 



/l-4a" / l-2a" Q Q; 

/ ^ ff (a,/3)ff(/3, 7 )ff(7,«")rf7^, YTQa~ < < 1 + 4^ 

1+2^ 

f 8 (E ,E) = 0, <a "<—— (3.7) 
All these integrals can be solved analytically. 



3.6 High Energy Gamow-Teller states 
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In this same paper, Jin gives the parameters A4 and Bi as a function of the 
energy. 

A l = -0.937680 + 0.203032 x lO 1 ^ ' 2 - 0.892372E " a4 



M = -0.110826 x 10- + xl<r> 



0.185895 

E 1 
^0 



= 0.721355 x 10" 



A 4 = -0.775913 x 10 1 + 0.204322 x lO 2 ^ 1 - 0.124189 x 10 2 £°' 2 



A 5 = -0.343254 x 10~ 2 + 0.213600 x 10 _1 £ - 0.882892 x 10 -3 ^ 

B 2 = -0.156849 + 0.965065Eo , and 
B 3 = 0.443463 x 10" 1 + 0.1697733 x lO" 1 ^ 

As a first attempt to fit the calculated background we used the 1407 keV peak. 
This peak was chosen because its Compton edge can be easily distinguished in the 
54 Ni and 54 Co spectra. The integrals in (3.7) and (3.8), were solved analytically with 
MAPLE10* m and its solutions with the full energy peak, single Compton scattering 
and the flat continuum are presented in Figure 13.151 

In Figure 13". 161 is shown the 54 Co total spectrum with the calculated background. 
A correction factor was applied to this calculated background in order to equal the 
amplitude of the full energy peak of 1407 keV with the measured spectrum. A good 
agreement between the theoretical background and the real data is presented, but 
no further analysis has been done up to the time of writing. 



3.6 High Energy Gamow-Teller states 

In Figure 13.171 we present the energy spectra obtained in the 54 Fe( 3 He,t) 54 Co 
charge exchange reaction [21] . In addition to the + — > + Fermi transition and the 
first + — > 1 + Gamow-Teller transition, other high energy states are populated in 
the reaction. Because in the (3 decay case we have a Q-value window up to ~8.8 
MeV we also expected to populate some of these states. Because of the expected 
drop in /3-branching due to the / factor and the drop of the 7-ray detection efficiency 
at high 7-energies, we expect to see the 3375 keV peak more easily than the other 
higher lying states. We therefore searched for the 3375 keV peak. 

In the Figure I3TT81 the 3375 keV region is plotted, but no 7 lines are distinguish- 
able from the background. 
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3.7 Final Ti/ 2 and Absolute Gamow- Teller calcu- 
lation 

With the three values obtained from the behaviour of the 937 keV 7 line the value 
of the 54 Ni halflife is calculated as the weighted mean H which value is Ti/ 2 = 114 ±3 
ms. With this value we can estimate the absolute Gamow- Teller strength, using 
equation 11.6.21 

B(F)(1-5 C ) ( F A 2 v ~ 



T 1/2 Ka F \ JF R 2 



-GT 




where T1/2 = 0.114(3) s is the measured halflife, B(F) = 2 is the \N — Z\ value, 
(1 — 5 C ) = 0.995 is the Coulomb correction factor, a F = 4905 ± 245 is the number of 
counts in the + — > + transition in the charge exchange reaction (see Figure 15. 171 ). 

fp = 22107 ± 746 and Ei= G r a ? T fi = ( 130 - 9 ± 8A ) x 1q6 > K = 6144 - 4 ± 16 where 
the / values are given by the /-function calculation, and A 2 = 1.603 ± 0.007. With 
these values we can calculate the R 2 value which is related to the B(GT) as is shown 
in Eq. (1.24). 

1 2 x 0.995 

x 



0.114 ± 0.003 (6144.4 ± 1.6) x (4905 ± 245) 

x 



o 1.603 ±0.007 , . f 

1.0844 ± 0.0908 x 10 8 + — x 130.9 ± 8.1) x 10 £ 

R 2 



R 2 = 8.56 ± 0.73 
fi(GT) (Tl/2=114±3 ) = 0.48 ± 0.08 . (3.8) 
The calculated absolute Gamow- Teller strength is B{GT) = 0.48 ± 0.08. 



2 This is a preliminar estimation of the halflife, and certainly is not the right way to do it because 
this three values are not calculated from independent measurements. 



3.7 Final Tu 2 and Absolute Gamow-Teller calculation 
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Plsstics 

MB Heidelberg MB Leuven 



Figure 3.3: Veto Condition: If a f3 particle is detected in one of the side plastic 
detectors, the corresponding 7 's are accepted if they happens in the opposed MINI- 
BALL cluster. If the (3 particle goes through the middle (3 detector, the 7 rays are 
accepted if they happens in the two most distant crystals of each MINIBALL cluster 



32 



Chapter 3 




Figure 3.4: Difference between the no-vetoed (up) and vetoed (down) laser on Ni, 
mass A =54 gated 7 spectra 
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Figure 3.7: Interest zone zoom of the A=54 laser on Ni and laser on Co spectra. 
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Figure 3.8: ID Curve during a macrocycle time (400 ms implantation and 600 ms 
decay) solving the differential equations system with MaplelO. 54 iVz is shown in red 
and 54 Co is shown in green. 
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Laser on Ni. 937 keV evolution on time. Gaussian fit with a linear background substraction 




Figure 3.10: Implantation- decay curve of 54 Ni. Each point represents the area of the 
fitted Gaussian, subtracting a linear background in the near region. 
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Laser on Ni. 937 keV evolution on time. Double Gaussian fit with a long linear background substraction 




Figure 3.12: Implantation- decay curve of 54 Ni. Each point represent the area of a 
Gaussian with fix width and centroid at 937 ke V, subtracting a long linear background 
taken from 800 to 1100 keV, in the region. 
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Laser On Ni. 937 keV evolution on time. Integration and background subtraction 




Figure 3.14: Implantation- decay curve of 54 Ni. Each point represents the area of the 
fitted two Gaussian subtracting a long linear background, taken from 800 to 1100 
keV, in the region. 
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Figure 3.15: Background estimated from Jin's semi empirical method. The pink line 
represents the single Compton scattering, the yellow line is the Gaussian peak, the 
green line is the flat continuum, the red line is the double Compton scattering and 
the blue line is the triple Compton scattering. The light blue line is the sum of all 
of them. 
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Figure 3.16: 54 Co Spectrum is shown in black and the semiempirical background 
calculation is shown in red. 
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Figure 3.17: Energy spectum of the Fe(^He,t) Co charge exchange reaction. 
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Figure 3.18: 5 ANi spectrum in the 3375 keV vicinity. 



Chapter 4 

Conclusions and Future 



The change exchange reactions and the /3-decay process are complementary ways 
to study the Gamow- Teller strength. The advantage of the charge exchange reac- 
tions are that they are not limited by the Q-value as the /3-decay is. But the big 
advantage of the /3-decay is that the measurement of the Gamow- Teller strength to 
one particular state is sufficient to normalize the Gamow- Teller strength obtained 
by the charge exchange reactions. 

In this report a /3-decay experiment was described and analyzed. The total halflife 
of the 54 Ni was obtained from the 7 spectra giving a tentative value of Ti/ 2 = 114 ±3 
ms. In order to improve and facilitate the halflife analysis, a semi-empirical model 
for the 7-ray response function of a Ge detector was used to obtain the background 
shape including multiple Compton scattering and the Compton edge. This phase of 
the analysis is still under development but we already obtained satisfactory results 
if we compare the shape of the estimated background with the real shape obtained 
in our experiment. 

In the near future we will start with the analysis of the (3s in order to obtain the 
branching ratio between the super-allowed Fermi and the first Gamow- Teller states. 
A new experiment to measure the decay of 54 Ni has been performed very recently 
at GSI (Darmstadt, Germany) and the on-line analysis showed very promising re- 
sults. The main advantage of this facility is the high rate of production of 54 Ni in 
the Fragment Separator (~100 54 Ni ions per second, in comparison to ~ 3600 total 
counts at LISOL). For that reason we hope to see the high energy Gamow- Teller 
state in the /3-decay, as well. 



Bibliography 



H. H. Heck-man & P. W. Starring, "Nuclear Physics and the Fundamental 
Particles". Holt, Rinehart and Winston, INC. (1963). 

R.R. Roy & B.P. Nigam, "Nuclear Physics: Theory and Experiment". 
John Wiley & Sons, INC. (1967). 

A. Bohr & B.R Mottelson, "Nuclear Structure. Vol 1: Single Particle 
Motion". World Scientific Publishing Co.Pte.Ltd. 2 nd Edtion (1998). 

D. D. Warner, M. A. B entity and P. Van Isacker, Nature Physics 2 311- 
318 (2006). 

K.S. Krane, "Introductory Nuclear Physics". John Wiley & Sons, INC. 
2 nd Edition (1988). 

T.N. Taddeucci et al, Nucl. Phys. A469 125-172 (1987). 

E. Hagberg et al, Phys. Rev. Lett. 73 396-399 (1994). 
Y. Fujita et al, Nucl. Phys A687 311c-320c (2001). 
R.G.T Zegers et al, Phys.Rev. C 74 024309 (2006). 
Y.Fujita et al, Phys.Rev. C 67 064312 (2003). 

CD. Goodman et al, Phys. Rev. Lett 44 1755 (1980). 
Y. Fujita. Journal of Physics: Conference Series 20 107-112 (2005). 
Y. Fujita. Journal of Physics: Conference Series 49 29-34 (2006). 
Y. Kudryavtsev et al, Nucl. Inst, and Meth. B 114 350-365 (1996). 
Y. Kudryavtsev et al, Nucl. Inst, and Meth. B 179 412-435 (2001). 
P. Van den Berg et al, Nucl. Inst, and Meth. B 126 194-197 (1997). 

I. Reusen et al, Phys.Rev. C 59 2416-2421 (1999) 

O. Ivanov. PhD thesis. Katholieke Universiteit Leuven. 



[19] /. Reusen. PhD thesis. Katholieke Universiteit Leuven, 1999. 

[20] Y.Jin, R.P. Gardner and Verghese, Nucl. Inst, and Meth. in Phyics 
Research A 242 416-426 (1986) 

[21] Y. Fujita et al. Private communication. 



54 Cobalt total spectrum | 




Energy keV 



CO 

H — 1 


5000 






o 


4000 


O 






3000 




2000 




1000 









> 

a* 

CD 



O 

ccs 



CD 

> 
_cd 

CD 



CO 

c 

CD 



1.2 

1 

0.8 
0.6 
0.4 
0.2 


3000 
2000 
1000 




(a) c 



in 
< 

CO 











CD 



50 



Crf He,t) 50 Mn 

E=140 MeV/nucleon 
9=0° 




OJ 
CD 
CO 

CO 



3 



W 



50 



E in Mn (MeV) 







I 


p-decay f-factor: bU Fe --> bU Mn 








Q =8.152 MeV : 








EC 


— i — i — 

(c) 


— i — i — h 
+ 

q 


— I 1 — 

+ 


H 1 1 1 1 1 1 I I I I I I I I 1 1 1 1 

p-decay intensity: 50 Fe --> 50 Mn 










< 


LO 
CD 


+ + 




CO 


O 


-r-" C\T 
-i- CD 
^ CO 








Cxi co 






^ l ■ 


'- ' * I ' ^ ^ I ' * 



50 



E in Mn (MeV) 



54 Nickel total spectrum j 




Energy keV 



54 Nickel total spectrum | 




400 600 800 1000 1200 



Energy keV 



54 Nickel total spectrum j 




i i i i I i i i i I i i i i I i i i i I i i i i I i i i i I i i i i I i i i i I i i_ 

700 750 800 850 900 950 1000 1050 1100 

Energy keV 





^_ — I |H.1.UB_. ■ L. — -LI- 1— . Lii WTL 



iliiiiiiiii iMmumilZL 



p± ML J U-l .L,L. U_ C^LL Jll L 

f I 

'kUi,, hi LiijuLi 




E 



UNIVERSITY" DE VALENCIA 



I I I 
1 J J. 
I I I 
I I I 



I I I I 

L 1 J J. 
I I I I 
I I I I 



_ L 



I I I 
1 J J 
I I 
I I 



_ L 



I I I 
1 J J. 
I I I 
I I I 



I I I I 
L 1 J J. 
I I I I 
I I I I 



I I I I 
LI J J. 
I I I I 
I I 



I I I I 
_L1 J J. 



I I 



J J 
I I 
I I 



I I I 
J__l J 
III 
I I I 

Tin 
i i i 

Til' 

I I I 

"I — I" — I " 

I I I 

-I — I- — I ■ 

I I I 
f 4H' 

I I I 
4-4-1 

I I I 

J I _l . 

I I I 

1 J J. 

I I I 

I I I 



I I 
J_ J 
I I 
I I 
Tl 
I I 



I I 
I I 



I I I 

I J J 

III 

I I I 

T 1 

\y( i 


i i i 

r -t -i 

i i i 

-i — i- — i - 

i i i 

H — I — I ■ 
I I I 
4-4-1 
I I I 
J. 4 — I . 
I I I 
1 J J. 
I I I 
I I I 



I I 
I I 



!_ j_ 
i i 
i i 
r t 
i i 
r t 
i 

r 
i 

+ 
i 

4- 
i 

4- 
I 

.. 

I 

1 
I 
I 



I 



I I 
I I 



- r 



- h 



- I- 



- U 



_ l_ 



_ L 
I 
I 



_ 

I 

1 
I 

-t 
I 

4 
I 

4 

I 

_ 

I I 
I I 



- r 



i 

_ i_ 
i 

_ L 
I 
I 



I I 
I I 
Tl 
I I 
m 
i i i 
r t -i 
i i 

-\ 
i i 

4 -4 

i i i 

4-4-1 

i i i 

J. 4 _l 

I I I 
1 J J 
I I I 
I I I 



- r 



- h 



- 1- 



- u 



_ l_ 



_ L 



- r 



- h 



- I- 



- U 



_ l_ 



_ L 
I 
I 



I I 

I I 
1 T 

I I 
11' 

I I 
t T 

I I 
-t -\- 

I I 
4 —I ■ 

I I 
4 -I. 

I I 
4 

I 

J 
I 
I 



_l_ 



_l_ 



I 

- r 



- h 



- I- 



- U 



_ l_ 



_ L 
I 
I 



-|- 



-I- 



-I- 



_l_ 
I 



J J. 
I I 
I I 



I I 
I I 

r t 
i i 
r t 
i 

r 
i 

+ 
i 

+ 
i 

4- 

I 

I 

1 
I 
I 



I 

- r 
i 

- h 
I 

- h 
I 

- U 
I 

_ L 
I 

_ L 
I 
I 



I I 

I I 
"TT 

I I 
11' 

I I 
~t ~ l 1 

I I 
-\ —\ ■ 

I I 
4 H ■ 

I I 
4-1. 

I I 
4 _l. 

I I 
J J. 

I I 

I I 



I I I 

I I J 
I I I 

I I I 

Tin 
i i i 

T 1 
I I 

r 
i 

+■ 
i /i i 

4 -I 
I I 
4 -I 
I I I 

j I 



I /I I I 

1 1 J 
I I I 

I I I I 

r t i - i 
i i i i 
rni 
i i i i 
rni 
i i i i 
i — -i— — i- — i 
i i i i 

1-4-4—1 
I I I I 
I-4-4-I 
I I I I 
I i i 



I I I 

J__l J 

I I I 

I I I 

Tin 
i i i 

Til 
I I I 

r -r -i 

i i i 

-i— — i- — i 

i i i 

4-4—1 
I I I 
4-4-1 
I I I 
I i i 



I I I 

1 J J 
III 

I I I 

Tin 
i i i 

Til' 
I I I 

r t -i 
i i i 
4- -n 
i i i 

4-4—1 
I I I 
4-4-1 
I I I 
l l l 



I I 

!_ j_ 
i i 
i i 
r t 
i i 
r t 
i 

r 
i 

+■ 
i 

4- 
I 

4- 
I 
I 



I I 

!_ j_ 
i i 
i i 
r t 
i i 
r t 
i 

r 
i 

+■ 
i 

4- 
I 

4- 
I 



I I 
I I 



I I 
J J 

i r 
i i 

in' 
i i 
~r ~r 
i i 
i -r 
i i 
-t —\ ■ 
i i 

4 H 
I I 

4-1. 
I I 

J I 



— r ~t /i- 



- r 



- h 



- h 



- u 



- r 



- h 



- 1- 



- u 



i 

- r 



- h 



- h 



- u 



—i — 



i 

- r 



- h 



- 1- 



- u 



-i- 



i i 
i i 
r t 
i i 
r t 
i 

r 
i 

+ 
i 

4- 

i 

4- 

i 



i 

- r 



- h 



- h 



- u 



i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i 
I I I I | I I I I | I I I I | I I I I | I I I I | I I I I | I I 







100 



200 



300 



400 



500 



600 



700 



Implantation Curve 



937 keV ID curve. Double Gaussian Fit. 
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